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Abstract. The rcnormalized trajectory in the muhi- dimensional coupling parameter space 
of the two-dimensional 0(3) non-linear sigma model is determined numerically under 
5-function block spin transformations using two different Monte Carlo renormalization group 
techniques. The rcnormalized trajectory is compared with the straight line of the fixed point 
trajectory (fixed point action) which leaves the asymptotically free ultraviolet fixed point 
of the critical surface in the orthogonal direction. Our results show that the renormalized 
trajectory breciks away from the fixed point trajectory in a range of the correlation length 
around § » 3-7, flowing into the high temperature fixed point at g = 0. The analytic large N 
calculation of the renormalized trajectory is also presented in the coupling parameter space 
of the most general bilinear HamUtonians. The renormalized trajectory in the large JV ap- 
proximation exhibits a similar shape as in the N = 3 case, with the sharp break occurring at 
a smaller correlation length of § « 2-3. 
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1. Introduction 

One of the most important issues in lattice field theory is the removal of cutoff effects in 
physical quantities. Two different approaches have recently received considerable attention, 
with the promise of efficient reductions in cutoff contaminations. The first approach, which 
combines the original Symanzik improvement program [1, 2] with tadpole improvement in 
the lattice coupling constant [3], has been successfully tested in the spectroscopy of heavy 
quark-antiquark bound states [4]. The second approach [5], which attempts to approximate 
the renormalized trajectory (RT) within the context of Wilson's renormalization group pro- 
gram [6], shows considerable promises for lattice QCD applications [7, 8]. In this work, we 
determine the RT in the two-dimensional 0(3) non-linear sigma model and compare it with the 
fixed point approximation, as proposed in the second approach [5]. The precise connection 
between Symanzik's improvement program and a truncated approximation to the asymptoti- 
cally free fixed point in Wilson's renormalization group approach remains an interesting and 
yet imresolved issue. 

It has been known for a long time that, in contrast with the standard lattice action, each 
point on the RT defines a perfect lattice action [5] which is free of cutoff effects at any finite 
correlation length [6]. A sketch of the RT in the phase diagram of asymptotically free field 
theories is shown in Fig. 1 as plotted in the infinite dimensional space of couplings Ki. The 
inverse coupling 1 /Ki of the standard lattice action is singled out to label the horizontal axis. 
The critical manifold is given in Fig. 1 by the Ki = oo plane and the asymptotically free fixed 
point on the critical surface is designated by UVFP. The continuum theory at finite lattice cor- 
relation lengths is defined by the RT which flows along the unstable direction from the UVFP 
to the high temperature fixed point (HTFP) with vanishing correlation length atKi = [6]. The 
determination of the perfect lattice action has not been considered feasible in the past, mainly 
because a practical approximation to the RT in the infinite dimensional coupling parameter 
space was thought to require the uncontrolled truncation of a large number of terms in the 
blocked Hamiltonian. 

Considerable progress has been made recently by Hasenfratz and Niedermayer [5] who re- 
alized that the fixed point lattice action of the UVFP can be determined from a classical saddle 
point problem in asymptotically free field theories. The fixed point trajectory (FPT), which is 
defined as the straight line that originates from the UVFP and is perpendicular to the critical 
surface, has been suggested to be a good approximation to the RT at sufficiently large corre- 
lation lengths [5] (dashed line for if i < oo in Fig. 1.) It was also demonstrated that the fixed 
point action (FPA) can be rendered very short-ranged by optimizing the block-spin renormal- 
ization group transformations with the expectation that an approximate FPT which is based 
on a truncated FPA, with the very small long range couplings neglected, will exhibit almost no 
cut-off dependence at large and moderate correlation lengths. 

The performance of the approximate FPT has first been tested by Hasenfratz and Nieder- 
mayer in the 2d non-linear 0(3) sigma model which is known to be asymptotically free [5]. 
They carried out a pilot study using a truncated FPT with 24 different couplings. Several tests 
showed that the residual cutoff effects were not visible even down to a correlation length of 
three suggesting that truncation effects are negligible. This indirectly implies that the FPT 
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Figure 1 . Phase diagram of asymptotically free field theories. The inverse cou- 
pling l/Ki of the standard action, or equivalently the inverse correlation length, la- 
bels the horizontal axis, and the ratios Ki I Ki label the remaining (infinitely many) 
axes. The abbreviations RT, FPT, UVFP, and HTTP stand for renormalized trajec- 
tory, fixed point trajectory, ultraviolet fixed point, and high temperature fixed point, 
respectively. 

runs close to the RT in an extended range of lattice correlation lengths. 

For a better understanding of the FPT approximation, it is important to determine the cross- 
over region in the correlation length where the RT eventually has to break away from the FPT 
to flow into the HTTP as depicted in Fig. 1. In this work, we determine numerically the posi- 
tion of the RT in a finite dimensional subspace of the infinite dimensional coupling parameter 
space and compare it with the FPT. A more detailed account of our investigation will appear 
elsewhere [9]. 



2. The non-linear 0(3) sigma model and the fixed point action 

The path integral of the 0(N) non-linear sigma model in the continuvim is defined by 

Z = |D[<^(x)]exp[-;SJfcont] , (2.1) 

^cont = ^^d^xd„(t)(x)d^(t>(x) , (2.2) 

where jS = l/T is the inverse temperature, 4>{x) denotes the iV-component scalar field of unit 
length, and J D[</)(x)] is the 0(N) invariant measure. Ifwe integrate out the momenta between 
the cutoff A and A/t, we find the two-loop renormalization group result 

""-'r^^^T^^OiT^), (2.3) 



dhiT 2tt (2tt)2 
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which implies that the model is asymptotically free for Af > 3 and the UVFP is located at T = 

(j8 = 00.) 

The lattice Hamiltonian is not unique. The only constraints come from the symmetries 

of the system and the fact that it ought to reduce in the classical continuum limit to the Hamil- 
tonian Jfcont of Eq. (2.2). With these constraints, the lattice Hamiltonian can be parametrized 
as 

= Jf2 +^4 +^6 + ■■■ , (2.4) 



^2 = -^XXp(r) (l-<^);,<^;,+r) , 



2 . 



(2.5) 



Ma = X c(x,3/,z,w) (l-(/);,</)y) (l-</)z<^T^), (2.6) 

x,y,z,w 

where J-Cz includes all terms that are bilinear in the field variable 4>x, and 5/4, . . . denote all 
the other terms which are quartic or higher order in the field variables. The smnmations in 
Eqs. (2.5) and (2.6) are over two-dimensional lattice vectors r, x, y, z, and w. The standard 
lattice action is given by the first two terms in Eq. (2.5) with n = (1,0), (0, 1), and Ki = 
Pp(ri). Physical quantities are measured in units of the lattice spacing a which is set to one 
for notational convenience. 

Let us now consider the following 5-function block spin transformation, 

exp[-li'Mict>';p',c',... )] = J P(<^', 0) exp [-/?:^(c/);p, c, . . . )] , (2.7) 

where a blocked lattice site x' is assigned to a 2 x 2 cell of sites x on the unblocked lattice. 
We put primes in Eqs. (2.7) and (2.8) to label the quantities on the blocked lattice. For details 
on how the FPA can be computed in asymptotically free field theories we refer the reader to 
Ref. [5]. The authors of Ref. [5] showed that the fixed point couplings p*{r) agree with the 
ones obtained earlier for the non-interacting model with Gaussian block spin transformation 
[10], 



,ipr 



ip) 



smHpi/2) 



y n 

ni,^=-oo (Pi + 2Trni)2 -I- (p2 + iTinzV {pi/2 + mn) 



(2.9) 
(2.10) 



We have determined the couplings for several lattice vectors r by evaluating Eqs. (2.9) and 
(2.10) numerically. The values of the fourteen largest couplings are given in Table 1. The in- 
teractions are very short-ranged, since the couplings decrease rapidly when the distance |r| 
between two interacting spins grows. The couplings of the type p{r),r = (n, 0) fall exponen- 
tially with n, |p((n, 0))| ~ exp[-1.45n] [5]. The authors of [5] also demonstrated that the 
couplings can be rendered even more short-ranged by using a soft block spin transformation. 
Quartic and higher order fixed point couplings can be also determined from their saddle point 
equation. 
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p*(r) 


n 


= (1,0) 


-3.42839 


rs 


= (2,1) 


+0.01060 


r2 


= (2,0) 


+0.74981 


rg 


= (5,0) 


-0.00909 


r3 


= (1,1) 


+0.32486 


no 


= (4,1) 


+0.00828 


n 


= (3,0) 


-0.16871 


ni 


= (3,2) 


+0.00502 




= (4,0) 


+0.03877 


ri2 


= (5,1) 


-0.00274 


re 


= (3,1) 


-0.01976 


ru 


= (6,0) 


+ 0.00217 


r? 


= (2,2) 


-0.01658 


■'l4 


= (4,2) 


-0.00100 



Table i . The fourteen largest couplings p*{r). 



3. The renormalized trajectory in the large N limit 

In the large JV limit, Hirsch and Shenker derived the following recursion relation for the blocked 
two-point function in momentum space [11], 

^ ZiGip'/2 + nl)sHp'/2 + nl) 
Ti ^q'.i G{q'/2 + nl) sHq' /2 + nl) ' 

where I is a vector with components equal to or 1, and designates the number of points 
on the lattice before blocking. The momentum space propagator 

G{p) = —-^ 2 (3.3) 

p(p)+m2 

is the saddle point solution of the large N expansion (in the large N hmit, jS remains finite 
when N is factored out from the bilinear Hamiltonian.) For a given )3 and p(p), the mass gap 
m is determined from the gap equation 

^lG(p) = l. (3.4) 
^ p 

Eq. (3.1) can now be iterated to determine the blocked propagators after repeated block-spin 
transformations. In the subspace of bilinear Harmltonians, we can obtain the couplings Pp{r), 
P'p'ir'), etc. from the inverse propagators by Fourier transformation, i.e., for r,r' =t= 0, 

We have iterated Eq. (3.1) numerically, using two different expressions for the unblocked 
propagator. The first choice originates from the bilinear Hamiltonian that includes the four 
largest couplings in Table 1, 

pip) = aiX2(l -cosp^) + «2X2(1 - cos(2p^)) 
^ A* 

+ a3[2(l-cos(pi + p2)) + 2(l-cos(pi-p2))] + a4S2(l-cos(3p^)), (3.6) 
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Figure 2. flow //nes of pirz) / p{ri) (a) and p{r^) I p{r\) (b) in the large N limit. 
The flow lines represented by the crosses start from an initial propagator that is 
given by Eq. (3.6). The lines represented by squares, full, and open circles start from 
the fixed point propagator (cf Eq. (2.10)). The fixed point ratios p* (rz) / p* (ri) and 
P* {1^3) / P* {1^1) are marked by horizontal dashed lines. 

with coupling parameters ai, az, «3 and 014. The second choice is the fixed point propaga- 
tor resulting from Eq. (2.10). In the first case, we can start the flow lines, which after several 
blocking steps will approach the RT, at different positions in the four-dimensional coupling 
parameter space by varying the five coupling parameters ai , . . . ,014 and jS. In the second case, 
the flow lines start on the FPT and it is interesting to see how the blocked couplings break away 
from the FPT after a few block spin transformations while tracing the exact RT. The correla- 
tion length § on the unblocked lattice is given by § = V(«i -1- Aaz + 2«3 + 9a4)/m^ for the 
first choice, and by 5 = -JTfm? for the second choice (using the fact that p(p) in Eq. (2.10) be- 
haves as p(p) ~ for small values.) The parameter has been determined numerically 
from the gap equation (3.4). 

As an example, we have displayed in Fig. 2 the flow lines of the ratios pirz)! p{ri) and 
p{r3)/p{ri) as a fvmction of ln(l/5). Since the correlation length is reduced by a factor of 
two after every block spin transformation, two consecutive points on a renormalization flow 
line are always separated by an amount of In 2 on the horizontal axis. The flow patterns of 
three-link and other couplings look very similar and we refer the reader for more details to 
Ref. [9]. The UVFP is located in this plot at ln(l/§) = -00 and the HTFP at ln(l/§) = -1-00. 
The ratios p* (rz) / p* (ri) and p* (rs) / p* (ri) of the fixed point couplings are represented by 
horizontal dashed lines. The flow lines which start from the four-coupling Hamiltonian are 
represented by crosses. Fig. 2 shows that those flow lines that start at large correlation lengths 
are attracted to the FPT within a few block spin steps indicating that the FPT is indeed a very 
good approximation to the RT at large correlation lengths. At small correlation lengths, how- 
ever, the flow lines approach a curve (RT) which is substantially different from the FPT. The RT 
will eventually flow into the HTFP where the two ratios p {rz ) / p (n ) and p (rs ) / p (n ) vanish, as 
can be shown by the high temperature expansion [9]. Fig. 2 shows that the sharp break from 
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the FPT occurs at a correlation length of 5 ~ 2-3. The three flow lines that start on the FPT at 
three different values of ^ are represented by squares, full, and open circles. These flow lines 
trace the RT quite accvirately, since they start at a correlation length larger than three. 



4. Numerical simulation of the 0(3) model 

The renormalization group flows in the previous section were restricted to the subspace of bi- 
linear Hamiltonians in the large N limit. In contrast, the Monte Carlo renormahzation group 
method wfll allow us to study the exact flows of selected couplings in the 0(3) model. Al- 
though interaction terms which do not fit on the finite lattice are truncated, their effects are 
expected to be negligible. The only practical limitation of the method is the signal to noise 
ratio in measming very small blocked couplings. 

The 1-cluster algorithm [12] has been used to generate imblocked spin configurations on 
a 256 X 256 lattice using a Hamiltonian that includes the first fovu' couplings of Table 1. The 
spin configurations have been blocked four times down to a lattice of size 16x16 using the 
block spin transformation of Eq. (2.7). We have implemented two different methods to infer 
the couplings from the blocked spin configurations, a microcanonical demon method [13], and 
a canonical one [14]. The microcanonical method has been used earlier to determine blocked 
couplings in the two-dimensional 0(3) sigma model [15, 16], but for a block spin transfor- 
mation that differs somewhat from the one we are using in this paper. The microcanonical 
method performs very well from a computational viewpoint, with the disadvantage that the 
blocked couplings have a finite volimie bias [13]. In contrast, the canonical method is ex- 
act, but suffers from larger autocorrelation effects which make the method very cumbersome 
from a computational point of view. For more details about the simulation and a comparison 
of the two techniques, we refer the reader to Ref. [9]. 

To iUustrate our 0(3) results, we have displayed in Fig. 3 the renormahzation group flows 
of the ratios p(r2)lp(ri) and p{r^)/piri) as a function of ln(l/§). The correlation length on 
the unblocked lattice has been determined from the propagator in coordinate space by using 
the improved estimator technique of ref. [17]. Most of the renormalization group flow lines 
exhibited in Fig. 3 have been obtained with the microcanonical method (crosses.) We have also 
included renormalization group flow lines (circles) which were generated by the canonical de- 
mon method. For one of the canonical flow lines, the initial couplings were chosen to coincide 
with the starting point of a microcanonical flow line. Fig. 3 shows that the two flow lines with 
the same initial couplings coincide within the error bars, and the flow pattern is very similar 
to the one seen in Fig. 2. The only difference is that the sharp break from the FPT occurs now 
at a somewhat larger correlation length. The grey lines in Fig. 3 mark the approximate posi- 
tion of the RT, as determined from the flow lines. For comparison, we have also included the 
RT in the large JV limit (heavy fines) (cf. Fig. 2.) 



8 




ln( l/f ) 



FIGURE3. Flow lines of p{r2)/p(ri) and pir^) / p{ri) forN = 3. All flow lines start 
from the four-coupling Hamiltonian. The lines which are represented by crosses 
(circles) have been generated with the microcanonical (canonical) method, respec- 
tively. The grey lines indicate the approximate position of the RT, while the solid 
lines depict the projections of the large N RT onto the two planes. The fixed point 
ratios p* (r-i) I p* (ti) and p* {r^) / p* (ti) are marked by horizontal dashed lines. Er- 
ror bars are omitted when they are smaller than the symbol sizes. 

5. Conclusion 

We have demonstrated, by explicitly computing the RT in the 0(3) model, that the FPT pro- 
vides a good approximation to the RT in the region of the coupling parameter space where the 
correlation length is large. A significant break occurs, however, in a range of the correlation 
length around 5 ~ 3-7 where the RT sharply departs from the FPT and flows into the HTTP. 
The RT in the large N limit exhibits a similar behavior with the sharp break shifted to some- 
what smaller correlation length. Although we have only shown the projection of the exact RT 
in the K1-K2 and K1-K3 planes, several other projections of the flows were determined as well 
[9]. 

Both ends of the RT can be determined without the utilization of the Monte Carlo renormal- 
ization group. The RT in the large correlation length regime, near the UVFP, is well approx- 
imated by the FPT, while the RT in the regime of small correlation lengths, near the HTFP, 
can be obtained from a high temperature expansion. It is expected that both approximations 
will break down in the crossover region where the Monte Carlo renormalization group tech- 
nique may remain the only useful tool. It is not clear whether this regime will be important 
in practical applications. 

The studies we presented here will be useful to extend to non-Abelian lattice gauge theories 
in four dimensions [7]. 

Note added. After the completion of this work, M. Okawa pointed out some earlier results [18] 
where a similar large N approximation was applied to block spin transformations that differ 
somewhat from the ones used in our paper. 
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